Does quasi-long-range order in the two-dimensional XY model really survive weak 

random phase fluctuations? 
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Effective theories for random critical points are usually non-unitary, and thus may contain rel- 
evant operators with negative scaling dimensions. To study the consequences of the existence of 
negative dimensional operators, we consider the random-bond XY model. It has been argued that 
the XY model on a square lattice, when weakly perturbed by random phases, has a quasi-long-range 
ordered phase (the random spin wave phase) at sufficiently low temperatures. We show that in- 
finitely many relevant perturbations to the proposed critical action for the random spin wave phase 
were omitted in all previous treatments. The physical origin of these perturbations is intimately 
related to the existence of broadly distributed correlation functions. We find that those relevant 
perturbations do enter the Renormalization Group equations, and affect critical behavior. This 
raises the possibility that the random XY model has no quasi-long-range ordered phase and no 
Kosterlitz-Thouless (KT) phase transition. 

PACS numbers; 05.20.-y, 05.70.Jk, 75.50.Lk 
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I. INTRODUCTION 



The theory of the Renormalization Group (RG) and, in 
particular, the concept of relevant and irrelevant opera- 
tors provides a general and deep understanding of critical 
points in clean systems. The physical picture at the heart 
of the RG theory is of sufficient generality to have been 
applied to the study of critical points in disordered sys- 
tems such as random magnetic systems or the problem 
of Anderson localization . 

However, effective theories for critical points induced 
by disorder (in short, random critical points) are usually 
non-unitary, i.e., they contain operators whose scaling 
dimensions need not be bounded from below. Thus, it is 
possible for a random critical point to be endowed with 
operators whose scaling dimensions are negative. This 
possibility is closely related to the fact that observables 
can be very broadly distributed in a critical phenomena 
induced by disorder. 

A paradigm of this situation is the problem of a "rel- 
ativistic" particle moving in two spatial dimensions in 
the background of a static but random vector potential, 
in short the random-Dirac-fermion problem [Q-|lO| . This 
system was found to have a line of critical points such 
that each critical point contains operators with negative 
scaling dimensions that carry trivial quantum numbers 
associated to the symmetries in the problem In 
other words, those operators are relevant and may ap- 
pear in the effective theory. Now the question is, do the 
random-Dirac-fermion critical points really exist or are 
they destroyed by the relevant operators? To answer this 
question, the following issue needs to be addressed. Do 
these unusual relevant operators affect the critical points 
in the conventional way, can we use the standard RG ar- 
guments to study the stability of random critical points? 



One concern is that negative dimensional operators might 
have some special properties that would prevent the use 
of standard RG arguments. For example, whereas the 
identity operator has zero scaling dimension, it does not 
destabilize a critical point since it cannot affect the scal- 
ing of any correlation functions. The question we want to 
address in this paper is what happens to a critical point 
and the corresponding correlation functions when rele- 
vant operators with negative scaling dimensions appear 
in the effective theory. 

A second unusual property of the random-Dirac- 
fermion critical point relative to a generic critical point 
describing a clean system is that an infinite number of 
relevant operators appear simultaneously in the effective 
theory. The RG flow in the vicinity of the random-Dirac- 
fermion critical point must then involve infinitely many 
coupled equations in which case the issue of its stability 
becomes much more intricate. 

Since there is no general principle ruling out the ex- 
istence of an infinite number of relevant operators with 
negative scaling dimensions at a random critical point, 
it is imperative to reexamine the stability of random 
critical points with this possibility in mind. Usually, 
in the literature on random critical points, only the ef- 
fects of a finite number of perturbing operators are in- 
vestigated. The properties of the random-Dirac-fermion 
critical points suggest, however, that it is very impor- 
tant to study the scaling properties of "complex" op- 
erators. At a random-Dirac-fermion critical point, the 
scaling dimensions of those complex operators are of the 
form h = hon — gn^, where n is an integer that character- 
izes the complexity of the operator and g characterizes 
the strength of the randomness. The trademark of these 
complex operators is that, whereas they become more ir- 
relevant the more complex they are (the larger n is) in the 
absence of disorder, this behavior is reversed with disor- 
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der. In fact, since n is unbounded from above, it is seen 
that the random-Dirac-fermion critical points {g > 0) 
always have an infinite number of negative dimensional 
operators. For weak randomness, only complex opera- 
tors with very large n can have negative dimensions. In 
early studies of random-Dirac-fermion critical points, the 
effects of such complex operators (strongly irrelevant in 
the absence of disorder but relevant in the presence of 
disorder) were not accounted for. 

In this paper we want to gain more insights into ran- 
dom critical points characterized by an infinite number 
of operators with negative scaling dimensions, elucidate 
the physical origin of such operators, and study how these 
operators manifest themselves in the standard RG treat- 
ment of the stability of a random critical point. To this 
end, we will reexamine the problem of the random bond 
two-dimensional XY model 



H 



XY 



cos ( 



0j -^y)] 



(1.1) 



Here, the angles < <j>i < 27t are defined on the sites i of 
a square lattice. The positive exchange couplings or 
spin stiffness and the real phases Aij are defined on all 
directed nearest-neighbor pairs of sites (ij) and are in- 
dependently distributed on the links (ij) with variances 
(mean values) gj (J) and gA (0), respectively. 

Rubinstein et al. have studied the random bond XY 
model at the Gaussian level with randomness in the phase 
only {gj = 0, gA > 0) |ll|] . The Gaussian approximation 
consists in replacing Eq. (1.1) by the contimmm limit 



(1.2a) 



where (p is vortex free whereas Q carries vortices. The 
probability distribution for the static random vector po- 



tential A^^ is also taken to be Gaussian, 



/I?[A,]exp(-^/d2y Al) 



(1.2b) 



where we adopt the summation convention over repeated 
indices from now on. Rubinstein et al. argue that the 
random bond XY model belongs to the same univer- 
sality class as the Gaussian model and they infer from 
the Gaussian model that, for any given but sufficiently 
small disorder strength gA, there exists a line of critical 
points ending at a KT-like transition. In other words, 
for fixed gA the KT phase diagram of the pure system 
[ p2t is preserved albeit with scaling exponents depending 
on gA |l|]. At the heart of their argument is an esti- 
mate for the disorder average of a two-point correlation 
function for an operator associated with vortices. The 
relevance/irrelevance of this vortex operator controls the 
KT-like transition. 



The manifold of random critical points found by Ru- 
binstein et al. is quite special. First, these are non-trivial 
random critical points since scaling exponents depend 
both on temperature and disorder strength. This prop- 
erty should be contrasted with that of a critical point 
for which the effect of disorder is fully accounted for 
by irrelevant random perturbations as happens in the 
two-dimensional Ising model with weak bond random- 
ness p^ . Second, each random critical point is exactly 
soluble, all local operators can be listed and their scaling 
dimensions can be calculated. Third, there is an infinity 
of operators associated to vortices that carry negative 
scaling dimensions. We want to use the Gaussian ap- 
proximation to the random bond XY model as a testing 
ground to gain some insights about special properties of 
random critical points associated to a spectrum of nega- 
tive scaling dimensions without lower bound. 

We close this introduction by pointing out that, be- 
sides the relevance of the two-dimensional random XY 
model to magnetic systems with random Dzyaloshinkii- 
Moriya interactions |ll[], crystal systems on disordered 
substrates fl^ , arrays of Josephson junctions with posi- 
tional disorder ||l^, and vortex glasses ||T^, the random 
XY model is also closely related to spectral properties of 
two-dimensional Dirac Hamiltonians with random vec- 
tor potential and random mass. In turn, random Dirac 
fermions in two dimensions can be connected p7[ ] to sta- 
tistical problems such as the random flux-line model in 



18[ and driven ran- 
lat a better under- 



the mixed phase of superconductors 
dom diffusion model ||l^ . We hope t 
standing of the random XY model might be useful to 
this class of problems. 



II. RESULTS 

It is easy to show that, if vortices are not allowed 
(the spin wave approximation) , the continuum model Eq. 
( |l.2[ ) is at a fixed point where the spin exp(iiy9) has alge- 
braic quasi-long-range correlations for any temperature 
and any disorder strength. This phase, the random spin 
wave phase, is an exactly soluble random critical point. 

In this paper we would like to reexamine the stabil- 
ity of this quasi-long-range ordered phase. The stability 
of this phase has been studied before. If one assumes 
that the ground state in the vortex sector is in the dipole 
phase and if one considers the binding and un-binding of 
the simplest vortices within the first non-trivial order of 
a fugacity expansion for the vortices jl^, one finds that 
the random critical point is stable only for a range of tem- 
peratures and disorder strengths bounded by the dashed 
line in Fig. 1. A more sophisticated approximation con- 
sists in treating a non-interacting gas of dipoles in the 
presence of disorder non-perturbatively in the vortex fu- 
gacity, in which case quasi-long-range order is present in 
the shaded area in Fig. 1 [Q. However, in this paper we 
find that: 
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1. For any temperature and any disorder strength 
of the random phases, randomness in the vortex 
fugacity generates an infinite number of relevant 
terms (most of them carry negative scahng dimen- 
sions) in the critical action of the Gaussian ap- 
proximation to the random bond XY model that 
describes the random spin wave phase. 



2. 



If the continuum model Eq. (1.2) is deduced from 
a random phase only (gj = 0, gA > 0) XY model 
on a lattice, then its vortex fugacity is necessarily 
random. 



3. Most importantly, the above relevant terms in the 
critical action describing the random spin wave 
phase do enter the perturbative RG equations to 
each order in the fugacity expansion, and com- 
pletely modify the RG flow at long distances. Thus, 
the relevant terms have the potential to cause an 
instability of the random spin wave phase. 

Let us call the simultaneous presence of relevant terms in 
a critical action and in the RG equations a perturbative 
instability. Thus, we may say that the random spin wave 
phase of the random bond XY model has a perturbative 
instability for any temperature and disorder strength, if 
the vortex fugacity is itself random. 

We would like to stress that when there are finitely 
many relevant terms, the perturbative instability implies 
the instability of the critical point. Since they enter the 
RG equations, the finite number of relevant perturbations 
either destroys the algebraic long-range correlations alto- 
gether or changes critical exponents. 

However, in our case, the perturbative instability cor- 
responds to infinitely many relevant terms appearing si- 
multaneously in the critical action and in the RG equa- 
tions. It is thus not completely clear to us what are the 
effects of such a perturbative instability, after one sums 
up an infinite number of contributions to the RG equa- 
tions from all relevant operators. By contrast, if all but 
a finite number of relevant operators can be switched off 
from the critical action, these remaining perturbations 
would completely alter the correlation functions at long 
distances. 

At the very least, the perturbative instability in the 
random bond XY model represents a new situation with 
regard to the issue of the stability of random critical 
points which must be addressed. The possibility of this 
new situation (i.e., the presence of an infinite number of 
relevant operators in the critical action) is closely related 
to the fact that the moments of a random variable, say 
the exponentiated energy of a dipole of vortices, need not 
be bounded if the random variable is sufficiently broadly 
distributed. Hence, random critical points need not be 
described by unitary field theories and scaling exponents 
need not be bounded from below. 

Korshunov was the first to argue that there might 
not be any quasi-long-range ordered phase in the random 
phase XY model. To this end, he introduced a sequence 



of local operators, Or;Ar(x), for the replicated random 
phase XY model labeled by the two index r and N. For 
given integer values of r and N < r, Or;jv(x) creates N 
vortices, each belonging to a different replica of the XY 
model, on site x. Here, r is the total number of replica. 
Korshunov found that for any strength of disorder, the 
scaling dimension of Or.Ni^) becomes negative upon an- 
alytical continuation to the r | limit, provided N is 
fixed and chosen sufficiently large. 

We recall that the replica approach identifies a given 
physical operator Q with a family or sequence of oper- 
ators Qr labeled by the total number of replicas r. In 
this paper, we identify the family of operators Or-N (la- 
beled by r) studied by Korshunov with a physical opera- 
tor. More precisely, we show how the family of operators 
Or-N can be induced in a physical way in the effective 
action describing the random spin wave phase. 

In the absence of vortices, the disorder average over 
the two-point function of a local operator has the form 
A\x — y|" in the random bond XY model. On the one 
hand, if the fugacity expansion is valid (as is the case 
in the clean XY model), then the exponent a and the 
coefficient A depend on the fugacity Y and have an an- 
alytic expansion around Y = 0. On the other hand, the 
breakdown of the fugacity expansion can have three con- 
sequences: 

1. A{Y) is not analytic around Y — 0, but a{Y) is. 

2. Both A{Y) and a{Y) are not analytic around Y = 
0. 

3. The critical behavior is completely changed or de- 
stroyed by the inclusion of vortices. 

By extending the Renormalization Group (RG) equa- 
tions to fourth order in the fugacity, we can show that 
the fugacity expansion breaks down according to scenario 
2, i.e., both the coefficient A{Y) and the scaling exponent 
a{Y) are not analytic functions of the vortex fugacity. 
Furthermore, if conventional RG arguments apply, one 
may then conclude that vortices change or destroy the 
critical line of the the random XY model for any tem- 
perature and any impurity strength according to scenario 
3. 

The paper is organized as follows. We first show in sec- 
tion HI that the two-point correlation function studied 



in [ [ll[ to construct the phase diagram has a very broad 
probability distribution, a fact expressed in the underly- 
ing critical theory describing the random spin wave phase 
by the presence of infinitely many negative dimensional 
operators. 

We then show that the random spin wave phase is per- 
turbatively unstable. Indeed, this instability manifests 
itself by the non-analyticity of the fugacity expansion for 
disorder averaged correlation functions due to the exis- 
tence of infinitely many negative dimensional operators 
in the underlying critical theory. A summary of our ar- 
guments is presented in section |^ with technicalities rel- 
egated to appendix H. 
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Another signature of the perturbative instabihty, as 
we show in section 0, is illustrated by the fact that any 
randomness in the spin stiffness of the random bond XY 
model induces infinitely many relevant perturbations to 
the critical theory describing the random spin wave phase 
in the Gaussian approximation. 

The relationship between the random bond XY model, 
the random bond Villain model, and the Gaussian ap- 
proximation is discussed in section VI. It is pointed out 



that the fugacity expansion already breaks down in the 
random bond Villain model. 

Our conclusions are followed up by two appendices. 
The first one discusses ground state properties and some 
probability distributions for correlation functions are cal- 
culated. The second one presents a detailed derivation 
of the perturbative RG analysis up to fourth order in 
the fugacity expansion. There it is also shown that there 
exists a one to one correspondence between correlation 
functions for the perturbations induced by a random fu- 
gacity within the replica approach of section and the 
contributions to the fugacity expansion of section 



III. BROADLY DISTRIBUTED CORRELATION 
FUNCTIONS 

In this section we construct the random spin wave 
phase and show that it is described by a manifold of ran- 
dom critical points. The central quantity of interest in 
the spin wave phase is the thermal correlation function 
for two spins. We characterize uniquely its probability 
distribution by all its moments. All moments are well- 
behaved. 

We then turn our attention to the vortex sector and, in 
particular, to the exponentiated energy of a pair of vor- 
tices of opposite charges (a dipole) in the background of 
the random vector potential. All moments of this expo- 
nentiated energy are calculated in Eq. ( 3.18 ) . This is the 
central result of this paper. In contrast to the spin wave 
sector, higher moments dominate the lower ones. This 



property is nothing but the signature of a log-normal 
distribution for the exponentiated dipole energy. Corre- 
spondingl y, th e energy of a dipole is Gaussian distributed 
in model (^]^) . Conversely, the critical theory describing 
the random spin wave phase must contain infinitely many 
operators with negative scaling dimensions that are asso- 
ciated with vortices in order to account for the Gaussian 
distribution of the dipole energy. 



A. Factorization into a spin wave and a vortex sector 

We begi n with the model in the continuum defined by 
Eq. (1.2a). We reiterate that is assumed Gaussian 
distributed with variance gA, i.e., that 



(3.1) 



P[A^](xexp(-— / d'^Al 



The justification for this assumption is that the precise 
shape of the probability distribution should leave critical 
properties unchanged as long as the probability distri- 
bution preserves the short-range nature of spatial cor- 
relations in the disorder. In particular, the tails of the 
prob ability distribution of non-compact support in Eq. 
( |3.l[ ) should not affect critical properties. 

To understand the role of the random vector poten- 
tial, it is convenient to decompose it into transverse and 
longitudinal components: 



(3.2) 



The advantage of this decomposition is that the partition 
function derived from Hq = J (fxHc in Eq- (l-2a) and 
the probability distribution in Eq. (3.1) both factorize 
since 



P[9, rj] cx exp 



29A 



(3.3) 
(3.4) 



Here, O is dual to O ^, and one must implement the 
constraint on the disorder that there be no zero modes: 



j (fx 6l(x) = 0, J (fx 'nix) = 0. 



(3.5) 



The am bigu ity in the choice of the integration constant 
in Eq. (3.2) is thus removed. 



B. The spin wave sector 

The consequences of the random vector potential on 
the spin wave sector are trivial. One can shift spin wave 
integration variables to 



if-rj. 



(3.6) 



In the absence of randomness in the spin stiffness, all 
correlation functions for exp(i(^)=exp(i(/3' -I- iry) can be 



^ Given 0, the dual is defined by 9;j0 
example, if = X/i*=i arctan 



■ In I 



(x-Xi)2 
(x-Xi)i 



9^0. For 
then = 



where the vorticities rrii are integer and 



lo is an arbitrary length scale. 
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calculated. In turn, the disorder average over 77 can be 
performed since the probability distribution for i] is that 
of a free scalar field in two dimensions. 
For example, 



|yi -y2|^ 
1 



|yi-y2| 



(3.7) 



Thermal averaging is denoted by angular brackets. Dis- 
order averaging is denoted by an overline and induces 
a quadratic dependency on the moment q for the scal- 
ing exponent. Thus, the impact of the quenched random 
vector potential on the spin wave sector is to drive the 
system to a new critical point for any strength of the 
disorder gA and for any reduced spin stiffness 



K := J/T. 



(3.8) 



The random vector potential is seen to destroy the long- 
range order at vanishing temperature by replacing it with 
quasi-long-range order. For all finite temperatures, the 
algebraic decays of the spin correlation functions are 
more pronounced due to the disorder. On the other 
hand, random spin stiffness remains irrelevant since it 
amounts to a random temperature (more formally, one 
verifies that, for any integer valued q > 0, [d^Lp'Y'^ is 
a strongly irrelevant operator everywhere along the spin 
wave critical line K > Q). Note that this argument is 
nothing but Harris criterion | p2[ in disguise. Finally, by 
choosing |yi — y2| sufficiently large, the two-point func- 
tion (exp[ic/j(yi)— i(^(y2)]) is seen to be a random variable 
with an arbitrarily small random component. This is not 
so on all counts in the vortex sector. 



ScB[Q:e]: = EY,{T^k-nkf 



(3.9) 



-ttK ^(mfe - nk){mi - ni) In 



of Eq. (3.2), provided gA is not too large and assuming 
the existence of a dipole phase of the CB gas at suffi- 
ciently low temperature and large reduced bare vortex 
core energy i? |fl 2^. Again, Q is induced by a neutral 
configuration of vortices with vorticities ruk € Z, whereas 
is induced by a neutral configuration of vortices with 
vorticities ni G R. It is sufficient to consider neutral con- 
figurations since the energy cost of creating net vorticity 
scales logarithmically with the system size L, Iq being an 
arbitrary length scale. 

The perturbative RG analysis in the CB gas represen- 
tation is usually summarized by the phase diagram of 
Fig. 1. The phase diagram is three dimensional with 
1/K ~ T/J the dimensionless temperature, gA measur- 
ing the disorder strength, and Yi = exp(— i?) the fugac- 
ity of charge one vortex. All points on the plane with 
vanishing fugacity are critical. This is the manifold of 
critical points describing the random spin wave phase. 
Critical points within the shaded area are argued to be 
stable |ll| , po| , i.e., Yi is irrelevant and thus decreases at 
long distances. Critical points outside the shaded area 
are unstable, i.e., Yi is relevant and thus grows at long 
distances. 
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C. The vortex sector 



Vortices in the XY model are described by the field 
Q. More precisely, the local density of vortices on the 
Euclidean plane is given by d^Q. Only the component 
0^9 of the random vector potential A^^ couples directly 
to the vortices described by 8. Whereas the field Q is 
induced by integer valued vortices, the quenched disor- 
der dfj,6 describes real valued vortices. Hence, the system 
tries to minimize the energy by screening the real valued 
quenched vortices with thermally excited integer valued 
vortices. However, by doing so, entropy is lost. The 
balance of energy and entropy could lead to a KT-like 
critical temperature separating a low temperature phase 
with positive free energy and a high temperature phase 
with negative free energy. 

In fact in the absence of randomness in the spin stiff- 
ness, the existence of a KT transition is suggested by a 
perturbative RG calculation in the Coulomb (CB) gas 
representation 




FIG. I. Proposed phase diagram for CB gas with quenched 
randomly fractionally charged vortices. 1/K is the reduced 
temperature, qa the variance of the Gaussian disorder dfj,9, 
Yi the charge one fugacity for thermal vortices. 

To go beyond these results (still keeping J non- 
random), we prefer the S ine- Gordon (SG) representation 
of the CB gas in Eq. ( |3.9| ) with ruk = ±1 only. By 
an expansion in powers of the "magnetic field" hi of the 
Boltzmann weight with Lagrangian 
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^sg[x,0] -.^hd^xf-^cosx+^xidle), (3.10) 



2t 



followed by an integration over we recover the grand 
canonical partitio n fu nction of the charge one CB gas 
derived from Eq. ( |3.9| ) provided one identifies 



N 



1=1 



(3.11) 



Notice that x couples to the disorder d^9 through a 
purely imaginary coupling, and that higher charges vor- 
tices mfc = ±2, • • • are easily incorporated with higher 
harmonic cos(fcx), fc € N. Hence, were it not for the 
"magnetic field" hi, x and Q would decouple after the 
shift of integration variable 



X 



X 



it 



(3.12) 



very much in the same way the spin waves decouple from 
the longitudinal realizations of the disorder. 

In the absence of disorder, one can establish the exis- 
tence of the KT transition by performing a perturbative 
RG analysis on the two-point function 



(-Pxi,X2)o-^ 



(3.13) 



I n sh ort, one first expands the right hand side of Eq. 
( 3.13 ) in powers of a very small fugacity hi/2t. Without 
a short distance cutoff a, all coefficients of the expansion 
in the fugacity are ill-defined. The arbitrariness in the 
choice of the short distance cutoff is used to derive RG 
equations obeyed by the fugacity and the reduced tem- 
perature. The RG equations are integrated to determine 
whether the initial assumption of a very small fugacity 
is consistent. The irrelevance, marginality, and relevance 
of the fugacity then determines the spin wave phase, KT 
transition, and disordered phase of the XY model, re- 
spectively. 

Rubinstein et al. [0 followed the same strategy in 
the presence of the quenched vector potential 8^9. More 
precisely, they performed a RG analysis of the fugacity 
expansion of two correlation functions: 



Gxi,X2 ^ ln(-Fxi,X2)5 



(^x„x2) 



/ 2?[x] e~^^^^^'^^ gix(xi)-ix(x2) 

2sGm ' 



(3.14) 
(3.15) 



to the first non-trivial order in the fugacity. Notice that it 
is necessary to include both Gxi,x2 and (i^xi,x2) to close 
the RG equations to the first non-trivial order in the fu- 
gacity. This is not surprising since taking the logarithm 
does not usually commute with averaging. 



The cru cial point of the fugacity expansions in Eqs. 
(|I|,|1|) is that every expansion coefficients depend on 
correlation functions calculated for vanishing "magnetic 
field" hi (fugacity hi/2t) such as 



(fyi ■ ■ ■ (i2y2„(ei[x(xi)-x(x2)+x(yi)+--x(y2„)])^^^ 



(3.16) 



on the one hand, but also such as 



(ei[x(xi)-x(x2)])^^^o/ j d^yidPy2€^[x(yi)-x(y2)]\ 

(3.17) 

on the other hand. For vanishing fugacity hi/2t, all aver- 
aged correlation fun ctio ns are algebraic and in particular 
[compare with Eq. (3.7)] 



|yi - y2| 
1 



|yi - yap'T-f^'iCi-ffAif?) 



(3.18) 



This is our most important result. We will make use of 
it in section IV to deduce that the fugacity expansion 
cannot yield analytic scaling exponents around vanish- 
ing fugacity. The remarkable property of Eq. ( |3.18 ) is 
that the scaling exponent becomes negative for any given 
temperature 1 / K and disorder strength ga as long as the 
moment q is sufficiently large. Hence, there must exist 
infinitely many local operators with negative scaling di- 
mensions in the effective theory describing the random 
spin wave phase. 



The correlation function in Eq. ( 3.18 ) is closely related 
to the strength of the interaction between two external 
charges in the CB gas. Indeed, as we show in appendix 



,ix(yi) 



ix(y2)\ 



/ii=0 



(3.19) 



Here, -ffi,2 is the bare (since hi — 0) energy of two vor- 
tices of opposite unit charges in the background of the 
vector potential d^O. The probability distribution of Hi 2 
is calculated in appendix ^ and shown to be a Gaussian 
distribution with a variance growing logarithmically with 
|yi — y2|- Hence, the random variable exp[(_fC/ J)i/i_2] 
has a log-normal distribution. For a fixed separation 
|yi ~ y2|j the random energy Hi 2 can take arbitrarily 
negative values as a consequence of our init ial assumption 
on the probability distribution in Eq. ( 3T ) . This fact ex- 



plains why the random variable (exp [ix(yi) ^ ix(y2)]) is 
unbounded from above. 

Correspondingly, the ratio of the q moment to the first 
one raised to the power q grows with |yi — y2| raised to 
the positive power -f 27rg^fC^g(g — 1), in sharp contrast 
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to the mo ment s of the logarithm of correlation functions 
[see Eq. ( AlO )] on the one hand, or to the moments of 
correlation functions in the spin wave sector on the other 
hand. 



We infer from Eq. (3.18) that a sufficien tly la rge mo- 
ment of the two-point function in Eq. ( 3.18 ) is not 
bounded from above for arbitrary large values of the sep- 
aration |yi — y2|. This property is a consequence of Hi,2 
being Gaussian distributed that can also b e understood 
as follows. On the second line of Eq. ( 3.18 ), the disorder 
average is dominated by realizations of the disorder with 
0{yi) very negative and 9{y2) very positive. Such config- 
urations are extremely rare for small separation |yi — y2|, 
since the cost {d^O)"^ will then be substantial, but they 
become more likely as |yi — y2| increases. 

The sign of the dependency of the scaling ex pone nts 
2-KKq (1 — qaKq) on the right hand side of Eq. ( p.l^ ) is 
thus the signature of broadly distributed random vari- 
ables. We will show below that operators with the same 
scaling dimensions appear in the replicated action if ran- 
dom spin stiffness is introduced. The same scaling ex- 
ponents have also been found by Korshunov to control 
the fugacity expansion of correlation functions such as 
CFxi.xa) to order 2q Q (see section IV and appendix 
y|). Korshunov concluded from this property of the fu- 
gacity expansion that any quasi-long-range ordered phase 
should be destroyed for arbitrary weak random vector po- 
tential dfj^O. An alternative conclusion, however, is that 
quasi-long-range order is characterized by scaling expo- 
nents that are non-analytic functions of the fugacity. We 
will come back to this scenario in section [V and in the 
appendices. 

The parabola 



9a 



















I TT K 



(3.20) 



(dashed line in Fig. 1) is obtained from the first moment 
g = 1 by requiring that the scaling exponent on the right 
hand side of Eq. ( 3.18| ) be "marginal", i.e., equals 4 [ pT| . 
It is argued to delimit the boundary between quasi-long- 
range order and disordered phase for 7r/4 < 1/K < 7r/2 

Instead of the reentrant phase transition for < 
1/K < 7r/4, the dilute vortex pair approximation pro- 
posed in 1^ suggests that the parabola should be re- 
placed by the dotted line in Fig. 1. The difficulties 
with the fugacity expansion are here bypassed altogether 
since it is possible to calculate the mean of G'xi.x2 non- 
perturbatively in the fugacity provided it is assumed 
that: (i) an insulating dipole phase exists and (ii) the 
interaction between dipoles can be neglected. Neverthe- 
less, it remains an open problem to show rigorously that 
at T = and for infinitesimal gA, the ground state con- 
figuration is in some quasi-long-range ordered phase con- 
sistent with assumptions (i) and {ii) (see appendix ^). 



IV. NON-ANALYTICITY OF THE FUGACITY 
EXPANSION 



We are now ready to describe the results obtained 
from the fugacity expansion of correlation functions in 
the SG model Eq. ( ^.1C| ). We will prove that the fu- 
gacity expansion is non-analytic. For the pure system, 
non-analyticity of the fugacity expansion is interpreted 
as the destruction of quasi-long-range order. In the pres- 
ence of a random vector potential, we cannot rule out the 
possibility that an exotic phase survives with quasi-long- 
range order characterized by scaling exponents which are 
non-analytic functions of the fugacity. However, even if 
a quasi-long-range ordered phase is present in the phase 
diagram of the SG partition function, we will show in 
section ^ that randomness in the spin stiffness induces 
infinitely many relevant perturbations to this critical be- 
havior. 

The mathematical reason for the breakdown of the 
fugacity expansion is that we are expanding a random 
function in powers of a random variable that takes val- 
ues outside the radius of convergence of the expansion. 
The physical reason for the breakdown is that the typ- 
ical ground state of the random phase XY model does 
not support long-range-order (the ferromagnetic state). 
For the KT transition to survive the presence of a weak 
random vector potential 8^,9, the typical ground state 
must contain a sufhciently large number of tightly bound 
pairs of vortices so as to destroy long-range order, but 
not sufficiently large so as to screen the bare logarithmic 
interactions of the vortices (see appendix ^). Since the 
vortex fugacity measures, to a first approximation, the 
density of vortices, a ground state with quasi-long-range 
order must imply the breakdown of a fugacity expan- 
sion around vanishing fugacity (the ferromagnetic ground 
state). 

Although the fugacity expansion is non-analytic, it is 
still useful to decide if non-analyticity reflects only that 
of scaling exponents of algebraic correlation functions or 
if it signals the breakdown of algebraic order. We have 
performed the fugacity expansion in the SG representa- 
tion to fourth order in the fugacity and could not, to this 
order, distinguish between an exotic algebraic phase from 
the complete destruction of quasi-long-range order. 

Lastly, the fugacity expansion is also instructive in that 
it allows to classify and understand the role played by the 
large numbers of local operators that can be constructed 
within the replica formalism. The close relationship be- 
tween the fugacity expansion and the replica formalism 
will be established below together with appendix |b| and 
section 

In principle, we would like to calculate the probabil- 
ity distribution of the two-point functions G'xi.x2 ^-nd 
(-Pxi.x2)- This is done in appendix ^ for vanishing fu- 
gacity. For finite fugacity, we are only able to calculate 
their moments perturbatively in the fugacity. We ex- 



pect Gxi.xs to be close to a Gaussian distribution since 



7 



it is already Gaussian distributed for vanishing fugacity. 



Hence, we will only calculate the mean of 



On the 



other hand, we will need all moments of (FX1.X2) since 
it is log-normal distributed for vanishing fugacity. Our 
goal is thus to calculate perturbatively in powers of the 
bare fugacity hi/2t 



Gi 



(4.1) 
(4.2) 



We restrict configurations of (thermal as well as 
quenched) vortices to neutral ones. This implies that 
only eve n powers of the bare fugacity hi/2t enter in Eqs. 
(4.l |4.2| ). The calculation to lowest order in the fugacity 
is summarized by Eq. (3. IS). 



A. Fugacity expansion up to second order 

Up to second order in the bare fugacity hi/2t, we find 
(see appendix ^) that the mean of the two-dimensional 
CB interaction between two external charges of opposite 
sign is 



G 



12 



27r 



L/a 



In 



X12 



(4.3) 



whereas 



X12 



-2'JTx{q) 



with the scaling exponents 



x{q):=K{q) 



K{q) 



L/a 



K'aW } ydi) J dyy 



(4.4) 



(4.5) 



3-27rif 



Here, we have introduced 



is renormalized multiplicatively, and X, K{q) axe renor- 
malized additively. The disorder strength qa is left un- 
changed in this scheme. We thus recover the well known 
RG equations [O 



dY, 



(la) 



= (2-7rX)y(i^i), 



di 



dK{q) 

dl 
dgA 



-A'K'^q ■ 



K{q) 



H 2 



y-2 



dl 



= 0. 



(4.8) 
(4.9) 

(4.10) 
(4.11) 



So far so good, the scaling exponent K t hat contro ls the 
algebraic decay of the two-point function {Fi2)hx=o gives 
us a criterion for the breakdown of the fugacity expan- 
sion. All scaling exponents x{q) are finite as long as 
4 — 2nK < up to second order in the bare fugacity 
hi/2t. This is not true anymore to fourth order in the 
bare fugacity. 

Let us stress a few important points. 

1. ^(1,1) appears in the RG equation of K. The latter 
controls the scaling dimensions of many operators. 

2. When 4 - 2TrK < 0, Y^i^) flows to zero. The RG 
flow of K [K{q)] only changes K [K{q)] by a fi- 
nite amount proportional to ^(i-i)- Thus, scaling 
exponents only receive corrections of order ^(i.i)- 
We may then say that the scaling exponents are 
analytic in ^(i:!) around y(i-i) = 0. 

3. When 4 — 2ttK > 0, ^(1,1) flows to infinity at long 
distances and K blows up. This was interpreted 
as the instability of the random spin wave critical 
point in Ref. ||ll|. At the very least, scaling expo- 
nents are not analytic in ^(i;!) around ^(i;!) = 
since the perturbative RG flow breaks down. 



Y2 



a? hi 
2t 



K:=K{1), Kiq):=Kq-K'g^q' 



(4.6) 
(4.7) 



The reason for which we label the dimensionless fugacity 
^(i;i) by the subindex (1; 1) will become clear to fourth 
order in the fugacity expansion. Suffices to say that to 
second order in the fugacity expansion, only one pair of 
vortices [of the type (1; 1)] renormalizes the moments of 
two-point functions. 

The short distance (dimensionless) cutoff in the ubiq- 
uitous integral on the right hand sides of Eqs. (4^,4.5) is 
arbitrary. By splitting the range of integration into two 
ranges [1, e' [ an d [e', oo[ with < Z ^ 1, it can be shown 
that Eqs. 



i) are form invariant provided Y( 



B. Fugacity expansion up to fourth and higher order 



We h ave performed the fugacity expansion of Eqs. 
( ^Tp^ ) to fourth order in the fugacity hi/2t. The fourth 
order calculation shows that, due to the disorder, the 
most singular coefficients in the fugacity expansions of 
6*12 and {Fi2)'^ are proportional to the integral 



lim 



dyy 



3-2itK(2} 



(4.12) 



Our derivation of the fugacity expansion of (-F12) is 
sketched in appendix It is also shown there that, 
in general, the most singular coefficient of the fugacity 
expansion is proportional to the integral 
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lim 



dyy 



3-2iTK(n) 



(4.13) 



(i.i;2) 
dl 



2 - ttK{2) 



Y, 



(i,i;2)- 



(4.16) 



to the order 2n. Hence, to the 2n-th order, the regime of 
analyticity of the fugacity expansion is dehmited by the 
curve [compare with Eq. (3.20)]: 



mi„[<,«(l//fl....«(l/A-):.ji;(l-2^ 

(4.14) 

in the plane of vanishing fugacity in the three dimensional 
couphng space depicted in Fig. 2. 

We must th erefore conclude that the fugacity expan- 
sions of Eqs. (4.1,4.2) are non-analytic for any infinites- 
imal value of the disorder strength qa-, since the region 
of analyticity in Fig. 2 shrinks to each new order and 
collapses to the segment < XjK < tt/2, gA — as 
n "f oo. Notice that this argument breaks down in the 
pure system where 



lim K{n) 

OA 10 



Kn. 



(4.15) 



This is nothing but the statement that cos(nx) is the 
more irrelevant the larger n is. To put it differently, it 
is necessary that vortices be more irrelevant the higher 
their charges for the fugacity expansion to be analytic. 




K 



FIG. 2. Boundaries g^2\^/K) g^2\^/K) in the plane 
of vanishing fugacity. The shaded area represents the regime 
of analyticity of the fugacity expansion to fourth order in the 
fugacity. 



The integral in (4.12) appears in the RG flow to fourth 
order in the fugacity expansion. Indeed, we find that: 

1. The RG equation of K contains a new fugacity 
^(i,i;2)- 

2. The RG equation for i^(i,i;2) is 



Therefore, when 4 — 2ttK{2) > 0, y(i,i:2) flows to infin- 
ity and the RG flow of K blows up, suggesting the new 
phase boundary 27ri4r(2) = 4, or, equivalently, g\ {1/K). 
We believe that a similar structure appears to the order 
2n in the fugacity expansion: 

1. The RG equation of K contains a new fugacity 
^(1, •■■,!;")■ 



2. The RG equation for y(i,...,i;„) is 



dl 



2-TrK{n) Y^^ 



(4.17) 



The quasi-long-range phase boundary is controlled by 
2TiK{2n) = 4, or, equivalently, g'^^\l/K). A sta- 
ble quasi-long-range ordered phase can only exist if 
2'KK{2n) > 4 for aU n € N. 

Based on the conventional RG analysis which essen- 
tially assumes that one can switch off all fugacities but 
to order 2n, we conclude that the quasi-long- 
range ordered phase is destroyed by any amount of dis- 
order. A more conservative conclusion that we can draw 
is that the perturbative RG flow must break down be- 
yond some critical order in the fugacity expansion that 
depends on the strength of the disorder. In any case, 
scaling exponents cannot be analytic functions of the bare 
vortex fugacity Yi . 



V. THE PERTURBATIVE INSTABILITY 



We have shown in section III that the critical theory 



describing the random spin wave phase must contain an 
inflnite number of operators with negative scaling dimen- 
sions. It was shown in section IV that these operators 
with negative scaling dimensions have dramatic conse- 
quences on the RG equations within the fugacity expan- 
sion. We now complete the proof of the existence of a 
perturbative instability by showing that the effective ac- 
tion from which disorder averaged correlation functions 
are built is necessarily perturbed by infinitely many op- 
erators with negative scaling dimensions. 

The SG representation makes it clear that random spin 
stiffness cannot be dismissed as irrelevant as was the case 
in the spin wave sector since it induces a random fugac- 
ity in addition to a random temperature. Random spin 
stiffness has two very different consequences from a sym- 
metry point of view. 

To see this, we fermionize the SG model. We use the 
(Euclidean) bosonization rules 



9 



on 



hi 



(5.1) 



-imiip^jj — > — — cosx, 



to relate bilinears of Grassm ann variables tp to the real 
scalar field x of Eq. (3.10). Here, 7^ are any two of the 
three Pauli matrices. Thus, the grand canonical parti- 
tion function of the CB gas derived from Eq. (3.9) is 
equivalent to first expanding the partition function with 
Thirring Lagrangian 



^Th:=^(i7M5M + i"^i)^-|(^7MV')' + (V57MV')(5^^) 



(5.2) 



in powers of the mass imi and then integrating over ^p, ip. 
The couplings g, imi are related to the reduced spin stiff- 
ness K by 



1 



hi 



(5.3) 



Randomness in the vector potential of the XY model 
thus enters the Thirring model through the coupling of 
the current := ip"f^ip of Dirac fermions to the trans- 
verse component 9^0 of the vector potential A^. Ran- 
domness in the XY spin stiffness enters the Thirring 
model through a random mass and through a random 
current-current interaction. Moreover, the random mass 
imi and random Thirring interaction g are generally not 
Gaussian distributed (even for a Gaussian distributed 
spin stiffness). 

The fermionization of the random CB gas shows that 
random spin stiffness has two very different consequences 
from a symmetry point of view. On the one hand it in- 
duces a random mass which breaks the chiral symme- 
try of the kinetic energy of the Dirac fermions. On the 
other hand it induces a random current-current interac- 
tion which preserves the chiral symmetry. Hence, any 
RG calculation should renormalize imi very differently 
from g. 

In fact, it is sufRcient to ignore randomness in g al- 
together for two reasons, provided we include random- 
ness in imi and 9^0. Firstly, randomness in g resembles 
non-Gaussian distributed vector potential 8^,9. It is then 
straightforward to show that non-Gaussian randomness 
in 8^9 is irrelevant by calculating the scaling dimension 
of ('07;jV')^' ~ i8fiX)^'^ on the critical plane of Fig. 1. 
Secondly, random imi is induced by a random fugacity 
of the CB gas (random magnetic field of the SG model) 
which will be seen to play a key role beyond the Villain 
approximation of the random phase XY model. Hence, 
it is meaningful to treat im, g and 9^0, as independent 
random fields and to assume that only randomness in im 
(or hi in SG model) and dfj,9 are present. 



To understand the interplay between a random fugac- 
ity Yi (magnetic field hi) and a random vector potential 
dfj^d, we use the replica formalism based on the identity 
Ina: = limrxo(a;'^ — 1)/?'. Indeed, within the replica ap- 
proach, disorder average can be performed directly on 
the replicated partition function. 

Thus, starting from the replicated Lagrangian 



a=l 



77(9^X0) - — cosxa + 7r-Xa[d^ 



2t 



r 



2tt' 



+Ja^id^e)\ , (5.4) 



0=1 



one can integrate over the disorder in hi (imi) or in 9^6'. 
In particular, integration over non-Gaussian hi (imi) 
generates r'', q £ N, interactions through terms such 
as 



^(cosXa)(x) 

.a=l 



(5.5) 



in the replicated Lagrangian. According to conventional 
RG arguments, the importance of these operators is mea- 
sured by their scaling dimensions at a given critical point 
labeled by 1/K and g A- One then verifies along the 
derivation of Eq. (3.18|) that, of all operators. 



{i^ii^iY (5.6) 



are the most and less relevant, respectively, on the critical 
plane mi = 0, l/if > 0, > 0: 



(Ol...l(x) Oi...i(0)) OC \^\~^-Kq(l-gAKq) ^ (5 7) 

(Oj(x) 0,(0)) oc |xp2^^'?'(i-9^^). (5.8) 

The crucial point we want to make is that for any value of 
gA, C'l-.-i becomes relevant as q is increased. Hence, con- 
ventional RG arguments predict that all critical points 
on the plane 1/K > 0, gA > are unstable with respect 
to random fugacity Yi (magnetic field hi). 

Next, we show that the existence of infinitely many rel- 
evant operators is not an artifact of an ill-defined replica 
limit. Indeed, there exists a very interesting line of crit- 
ical points at twice the bare KT transition temperature 
1/K = TT, gA > 0, Yi = in the phase diagram of Fig. 1. 
From the perspective of the Thirring model, it describes 
massless non-interacting Dirac fermions in the presence 
of a transverse random vector potential. All disorder av- 
eraged correlation functions of local operators were calcu- 
lated with the help of supersymmetric techniques appro- 
priate to non- interacting systems [0J|]. It was found in 
that each critical point has infinitely many primary 
ds and that these primary fields control the multi- 
fractal scaling of the random "wave function" exp[6'(x)]. 
These primary fields are precisely given by the family 



fie 



10 



Oi...i in the replica language and carry the scaling di- 
mension 2q[l — {gA/TT)q] appropriate to 1/K = n. 

If we accept the hypothesis that having infinitely (as 
opposed to finitely) many relevant perturbations results 
in the instability of the line l/A' = tt, > 0, we need 
not even rely on the replica approach to rule out a KT 
transition. Indeed, the shape of the phase digram of Fig. 
1 is preserved if we forbid thermal excitations of vorticity 
one (Yi = 0) but allow thermal excitations of vorticity 
two {Y2 > 0). More precisely, if we replace in the SG (Th) 
Lagrangian ^cosx (imi^ptlj) by ^ cos(2x) (im2('0'0)^)j 
then the parabola 5a(-^) = simply rescaled 

to -^(1— ^ X f-^). The inescapable conclusion is then 
that random fugacity for charge two vortices destroys this 
quasi-long-range ordered phase since it generates relevant 
interactions given by Eq. (5.7) under length rescaling. 



reached otherwise. Although the random phase can take 
all possible real values, the energy per link Cfj^ cannot 
take values outside the range (assumed compact for any 
finite temperature) 



< Cf,^ < 2 sup K,. 



(6.3) 



with probability one. 

The random Villain model consists in defining on the 
same lattice the random partition function pJ-pTl] 



V»=l / (u)'.j6Z 



(6.4) 



VI. THE VILLAIN APPROXIMATION 



Here, we have taken the spin stiffness to be selfaveraging, 
i.e.. 



In this section, we go back to the lattice to study in 
more details t he n ature of the Gaussian approximation 
made in Eq. (E2). We consider both the random bond 
XY and Villain models on the square lattice. We will see 
that the difficulties with the fugacity expansion are not 
associated with pathologies of the field theory at short 
distances (such as ill-defined operator product expan- 
sions) but are intrinsic to the fugacity expansion, i.e., 
are also present if the fugacity expansion is performed on 
the Villain model itself. We will also see that random- 
ness in the phase only {qa > 0, gj = 0) has different 
consequences in the XY and Villain models. 

We begin with the random partition function 



2L^ 



(6.1a) 



Cf^:= K,, [1 - cos(0, - c^j - A,,)] , (6.1b) 

for the random bond XY model on a square lattice made 
of sites. Directed links (two per site) on the square 
lattice are denoted by {ij). The phases Aij — —Aji are 
random (with short-range correlations for different links). 
However, they need not be restricted to < Aij < 2tt in 
spite of the periodicity of the cosine. Indeed, the proba- 
bility distribution for Aij need not be periodic with pe- 
riod 2tt. The reduced spin stiffness Jij > are also ran- 
dom (with short-range correlations for different links). A 
reasonable choice for the probability distribution of the 
random phases in the XY model is 



2L^ 



p{A,] ■■= n 



1 



(6.2) 



This choice is made without loss of generality provided 
any "small" departure from Eq. (^) does not pre- 
vent the system to flow to the fixed point it would have 



T \L 



(6.5) 



The periodicity under a shift of any by 2tt is preserved 
in the Villain action, but the non-linearity of the cosine 
has been removed in the Villain action. When referring 
to the random Villain model, we will have in mind the 
same probability distribution for the spin stiffness and 
for the random phase as for the XY model. 

The quant ity Qfj is not the counterpart to the link 



energy (a.lb) in the XY model since it is not periodic 



under a shift of Aij or (j^i by an integer multiple of 2n. It 
is, however, very closely related to the energy of a given 
configuration of vortices in the background of a random 
environment induced by bond randomness. 

By taking the random phase Aij of the Villain model 



to be Gaussian distributed according to Eq. (3.2), we 



immediately see that Qjj can take any arbitrary large 
value with a finite probability. Moreover, exp (Q^) is log- 
normal distributed very much in the same way as corre- 
lation functions for vortex operators are in the Gaussian 
approximation (see section III and appendix ^). This 
behavior should be contrasted with that of the Villain 
link energy 



:=-ln 



(6.6) 



which is indeed periodic under a change of Aij or </>,; by 
an integer multiple of 27r. It is crucial to realize that pe- 
riodicity of the Villain link energy £^ is broken to any 
finite order in a fugacity expansion, since the fugacity ex- 
pansion amounts to a truncation of the summation over 
lij €i Z. 

There is a noteworthy difference between the random 
bond XY and Villain models. If we take the spin stiffness 
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to be non-random but allow the relative phase between 
neighboring spins to be random, then we must assume 
that the SG model has both a random vector potential 
and a random fugacity if it is interpreted as the minimal 
model capturing the long distance properties of the XY 
model. To the contrary, the fugacity of the SG model is 
not random if it is derived from the Villain model with 
random phase but no randomness in the spin stiffness. 

To clarify these last two points, we need first to in- 
troduce some notation. We define the longitudinal and 
transversal components A| ■ and Aj^, , respectively, by 



where A' ■ is curl free, i.e.. 



0= curli A. 



(6.7) 



(6.8) 



--A 



-A 



-A 



■((i+x) ' ^^(i-|-x)(i+x-|-y) ' ^^(4+x+y)(i+y) 

and Ajj is divergence free, i.e., 
0= div, Aj-, 



-A 



II 

(j+y)i' 



A 



A 



i(z+x) ^(i — x)i 

Dual sites are labeled by 
1 



A 



i{i+y) 



—X H — y, 
2 2^' 



(6-9) 



(6.10) 



where the basis vectors x and y span the square lattice. 
For completeness, the gradient of a lattice scalar is de- 
fined by 



grad. 



M = x,y. 



(6.11) 



It is possible to rewrite the Villain partition function 
solely in terms of degrees of freedom defined on the dual 
lattice 




-ZyOC 



(ij> 



(6.12) 

We have dropped a multiplicative factor that depends on 
the (random) spin stiffness, and 



^i(i+A) 

(grad.(pi) 



qCB 



(6.13) 



2i(grad^¥>i) [At-(i+^^ + 2ttI^^, 



+/<) 



The random dual phases A^ are now curl free as are dual 
vortex degrees of freedom . 

Since Eq. ( 6.13| ) is quadra tic in the spin wave degrees 
of freedom (fi whereas Eq. ( |6.13 ) is linear in the vortex 



degrees of freedom l^, it is possible to decouple the spin 
wave sector from the vortex sector. Our final expression 
for the partition function in the vortex sector is [compare 
withEq. (U)] 



zcB^E n 

{mi} 1 j = l 

divi^ij =: — =="^i, 
V 27r 

diviAy =: \/27rni. 



-7r(mi-r!i)(mj-rij)_Dij[_ffkl] 



(6.14) 

(6.15) 
(6.16) 



Here, £'ij[iirki] are the random components of the dual 
lattice Green function in the background of random spin 
stiffness, namely the inverse of the quadratic form in the 
spin wave sector. 

Performing a dual transformation on the random bond 
Villain model thus offers two insights. First, given a 
Gaussian probability distribution for the random phases 
and no randomness in the spin stiffness, the probability 
distribution for the CB energy of a given configuration 
of vortices {mi} is Gaussian, and exponentiating this CB 
energy yields a log-normal distributed random variable. 
We thus conclude that the existence of log-normal dis- 
tributed correlation functions in the random spin wave 
phase is not an artifact of the continuum limit. 

Second, the diagonal components -Dii[_fi'ki] of the Green 
function define the vortex core energy. Hence, in the ab- 
sence of random spin stiffness, the core energy of the 
vortex sector is not random within the Villain model. 

We recover the XY partition function by replacing the 
right hand side of Eq. (6.4) by 



Mr 

^1 (2n)! 



E 



Aii 27rL; 



(6.17) 



Spin waves and vortices are coupled in the XY model 
with or without randomness. In an effective theory such 
as the SG model, non-linearities such as in Eq. (6.17) can 
be incorporated through a random fugacity to a first ap- 
proximation. In the clean limit, the irrelevance of higher 
vorticity charges justifies neglecting such an effect. How- 
ever, as we have shown in section ^ this is not anymore 
the case in the presence of random phases. 

A random core energy for vortices is always induced 
by a random spin stiffness. However, we have shown that 
the XY and Villain models differ in one very important 
aspect when only the relative phase between neighbor- 
ing spins is random. Indeed, due to non-linear effects the 
core energy of vortices is always random in the XY model 
with random phase whereas this is not the case for the 
Villain model. If randomness in the core energy does in- 
deed destroy the exotic quasi-long-range order proposed 
in [Q, we must then conclude that the Villain model 
with random phase only does not belong to the same 
universality class as the XY model with random phase 
only in a strong sense. It would be very interesting to 
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probe numerically our conjectured difference between the 
random XY and Villain models. 

We have also shown that the random XY model and 
any approximative treatment of the Villain model based 
on a perturbative expansion in the fugacity do not belong 
to the same universality class in a weak sense. We can 
illustrate this fact in a very suggestive way. We take the 
spin-spin correlation function (thermal average) in the 
random bond XY model for an arbitrary pair of sites i 
and j. Clearly, this is a random variable that takes values 
on a compact range, namely 

|(cos(/>,; cos(/)j)| < 1. (6.18) 

Consequently, all integer moments with p < q obey 



\{C0S (t>i COS (t)j)\P < 1. 



|(cOS0i COS0j)|9 < \{C0S COS <j)j)\P < 1- 



(6.19) 
(6.20) 



It is possible to verify that those inequalities are not sat- 
isfied to any finite order in a fugacity expansion of the 
random phase Villain model. Without loss of generality, 
we can take the continuum limit. We recall that the pre- 
scription to estimate the spin-spin correlation function 
(thermal average ) fr om the continuum theory of section 
[II is to identify [M 



.21) 



Zsc 



for some path joining Xj to x^. To any finite order in 
a fugacity expansion, the spin-spin correlation function 
(thermal average) is very broadly distributed (log-normal 
in the limit of in finite core energy), and thus violates the 
bounds of Eqs. (6.19 6.2C). We stress that this is a fail- 
ure of the fugacity expansion on the Villain model itself 
and not an artifact of the continuum limit ^. 

This situation is very reminiscent of that in the one- 
dimensional random bond Ising model. For a large but 
fixed separation, the Ising spin-spin correlation function 
(thermal average) is "close" to being log- normal |^ . The 
log-normal approximation describes exactly the mean 
and variance of the logarithm of the spin-spin correla- 
tion function (thermal average) in the Ising case but it 
neglects higher cumulants. By analogy, we might expect 
that the fugacity expansion on the Villain model cap- 
tures well the logarithm of the XY spin-spin correlation 



(thermal average). However, the log- normal approxima- 
tion in the random bond Ising model fails badly to de- 
scribe the tails of the Ising spin-spin correlation functions 
very much in the same way as the fugacity expansion on 
the Villain model dramatically overestimate tails for the 
probability distribution of the XY spin-spin correlation 
fimctions. 

So it is by now clear that the probability distribution 
of the random vector potential unduly favors rare events 
through its tails within the fugacity expansion. In one 
scenario, we must then anticipate partial loss of univer- 
sality for the fixed point probability distribution of corre- 
lation functions (thermal average) that are broadly dis- 
tributed, since their tails result from rare events. We are 
aware of several examples of this kind: directed polymers 
in a random medium |^], the metal- insulator transition 
p9|-^ , and quantum gravity jl^ . In the worst case sce- 
nario, the fugacity expansion on the Villain model is over- 
whelmed by rare events and loses complete predictability 
with regard to the phase diagram of the XY model. 



VII. CONCLUSIONS 



Concerned with the possibility that effective theories 
describing random critical points are often characterized 
by a spectrum of scaling exponents that is unbounded 
from below, we have studied the Gaussian approxima- 
tion to the random XY model on a square lattice within 
a perturbative RG framework. 

The Gaussian approximation to the random XY model 
predicts the existence of a manifold of random critical 
points describing a random spin wave phase. However, 
there are correlation functions in the random spin wave 
phase that are log-normal distributed. Correspondingly, 
there are infinitely many operators with negative scal- 
ing dimensions that are associated with vortices in the 
effective theory describing the random spin wave phase. 

We showed that all these operators with negative scal- 
ing dimensions contribute in a highly non-trivial way to 
the perturbative RG equations within a fugacity expan- 
sion around the random spin wave phase. The existence 
of infinitely many negative scaling dimensions thus man- 
ifests itself by the non-analyticity of the fugacity expan- 
sion for any given temperature and disorder strength. 



^An alternative way to stress this point is to note that the 
exact identity between the two- and four-points spin-spin cor- 
relation functions Pql 



^giW,-^,)) = j^gi(^,-^,)^|2 (5 22) 

that holds along the Nishimori line K = 1/<?a in the Villain 
model is always violated to any finite order in the fugacity 
expansion. 
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In this sense, the random spin wave phase is unstable, 
although we cannot preclude the possibility that a new 
phase with quasi-long-range order can be found for suf- 
ficiently low temperatures and weak disorder strength. 
The breakdown of the fugacity expansion is also associ- 
ated to a potential perturbative instability triggered by 
any random vortex core energy. 

We have shown that neither the breakdown of the fu- 
gacity expansion nor the perturbative instability are ar- 
tifacts of the continuum approximation that we used but 
would also be present in the random bond Villain model 
on the square lattice. Rather, they both reflect the ex- 
treme sensitivity of the fugacity expansion to the tails of 
the probability distribution that is chosen for the random 
bonds. 

The physical interpretation for the breakdown of the 
fugacity expansion is that the typical ground state is 
not ferromagnetic. In the best case scenario, the typi- 
cal ground state for weak disorder supports some quasi- 
long-range order that would persist for sufficiently low 
temperatures. However, to address the nature of the low 
temperature, weak disorder region of the phase diagram 
it is imperative to better characterize the typical ground 
state and to use a RG scheme that is non-perturbative 
in the vortex fugacity. 
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is ferromagnetic. In two dimensions, the long-range order 
of the ground state is first demoted to quasi-long-range 
order for any finite temperature below the KT transition 
temperature Tkt- In turn, topological excitations in the 
form of vortices wipe out the quasi-long-range order for 
T > Tkt- The crucial question we would like to address 
is: are typical ground states of 



{ij) 



(A2) 



ordered, quasi-long-range ordered, or disordered? 

What we mean by a typical ground state is the follow- 
ing. We assume that we know how to calculate the prob- 
ability to find a ground state with energy density e (en- 
ergy divided by total number of sites, i.e., < e < 4J). 
A typical ground state maximalizes this probability dis- 
tribution. Unfortunately, obtaining this probability dis- 
tribution is very difficult in view of: (i) the non-quadratic 
dependency of the energy on the disorder, and (n) of the 
need to minimalize the energy spectrum of 7J [{</), }; {^ij}] 
for a given realization of the disorder. The first difficulty 
can be disposed of in the Villain approximation (whereby 
it is assumed that the Villain approximation does not 
change the universality class), but we still must face the 
second difficulty. 

We focus on the vortex sector in the continuum ap- 
proximation of section [I] and show that, for extensively 
many realizations of the random vector potential 8^,9, 
the ground states are, for lack of a better description, 
complex configurations of vortices. We do this by gener- 
alizing an argument used in pc| ] to prove that the ran- 
dom potential d^d destroys the long-range order of the 
pure system at T = and for sufficiently strong disorder 
strength: gA > tt/S. We then modify the random phase 
XY model by restricting the possible realizations of the 
disorder to a more tractable subset. Within this subset 
we show that it is possible to decide whether the ground 
states support quasi- long-range order or not. 



1. Gaussian distribution of the energy of vortex 
configurations 



Let Q in Eq. (3.3) be given by the vortex configuration 



APPENDIX A: ZERO TEMPERATURE 
CONSIDERATIONS 



G(x) :— rrii In 



rrii G Z, 



(A3) 



In this appendix, we study some properties of the ran- 
dom phase XY model at zero temperature. Indeed, we 
recall that the prerequisite to the existence of a KT tran- 
sition in the disorder free XY model is the simple fact 
that the ground state of the Hamiltonian 



mm] :=E'^[i 



3(0i - 0j)] 



(Al) 



(which need not be neutral) and define its energy in the 
background of the random vector potential 9^0 to be 



(A4) 



where 
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^i,- -,n = -ttJ ^ m^m; In 

n 

SHi^...,„ = 27rjy^mfc6l(xfc). 

fc=i 



Xfe - X; 



^0 



(A5) 



The definition in Eq. (A4) is useful since the energy 
Hi^...^n is very closely related to the thermal average of 
the 2m — 71-point correlation function 

P — pix(xi) . . . ix(x™) -ix(xm+i) . . . -ix(x2m) 

Xi , ■ ■ ■,X2tti • C C C C , 

(A6) 

for vanishing fugacity hi/2t. Indeed, in that case 

(Fxi,.,x.„>hi=o = e-'^^i--^™. (A7) 

Notice that for a fixed realization of d^9, _ffi^....„ is 
bounded from below by — / d^x(9^6')^. On the other 

hand, since ^ J d^x{df^6y can take arbitrary large values, 
these are the rare events which take full advantage of the 
non-compactness of the Gaussian probability distribution 
for dfj,9, the probability distribution P(£'; Xi, • • • , x„) 
that Hi^...^n takes the value E is non- vanishing for all real 
values of E. A very special vortex configuration is the 
ferromagnetic configuration 9 = 0. This is the ground 
state of the pure system and it has a vanishing random 
energy. Clearly it need not be the ground state for a 
given realization of the disorder. 

What is the probability distribution of iJi _..._„? By 
definition it is given by 

P(£';xi,---,x„) 

JV[9] e'^/'''^ id.ofiy) ^(^- Jj-i,..,„) 

jv[e]e-^-f''"^'-^-'^"^^^ 



(A8) 



We can represent the delta function by an integral, in 
which case 



P(£^;xi, • • • ,x„) = 



27r 
dX 
2^ 



exp 



„iA(i5-ifi,....„) 
2(5Hi,...,„)2 



e 2 



(A9) 



Thus, the probability distribution for the energy of the 
vortex configuration 0(x) is Gaussian with mean i/i ... „ 
and variance 



((5i7i,...,„)^= -27rJ^gA ^ mfem/ln 



k,l=l 



Xfe - X/ 



lo 



2JgA Hi 



(AlO) 



We are now ready to calculate the (unnormalized) dis- 
tribution P{E; mi, • • • , m„) to find n vortices with vortic- 
ities mi, • • • ,m„ e Z anywhere on the Euclidean plane. 
It is defined by 



P{E;mi, - ■ ■ ,mn) : 
d^xi d^x„ 



P(£;;xi,---,x„), 



(All) 



where a is a microscopic cutoff and L is the system size. 
We are interested in the scaling of P{E; mi, • • • , nin) with 
the system size L. Notice that 

j dE P{E-mi,---,mn)=l^ . (A12) 
This scaling can be calculated in closed form for n = 1: 



exp 



P(S; m) 



In # 



VMgXHi 



Hi = m TrJln 



(A13) 
(A14) 



Here, the arbitrary length scale Iq is chosen to be L. 

In the absence of disorder, the ground state of Eq. ( A4 ) 
is the ferromagnetic state = 0. In the presence of dis- 
order, we c an c alculate the probability P(iJ;xi, • • • ,x„) 
from Eq. (A£) to find a vortex configuration with en- 
ergy E < 0. We say that it is energetically favorable to 
create a vortex configuration O different from O = if 
P{E < 0; Xi, • • • , x„) does not vanish. For example, to a 
good approximation, we find that 

f dE P{E;m)Ki P(0;m) 

J —00 

(A15) 



^^An'^J'^gAm'^ ln(i/c 



Equation ( ]A15 ) gives the number of sites on which it is 
energetically favorable to create a single vortex configu- 
ration. 

As proposed in ||2^], one can use estimate Eq. ( A15| ) 
to establish a criterion for the destruction of long-range 
order at T = of the pure system by the random vec- 
tor potential 9^6*. Long-range order is destroyed if the 
number of sites on which it is energetically favorable to 
create a single vortex configuration diverges in the ther- 
modynamic limit L 1 00, i.e., if the variance gA is larger 
than the critical value 



(A16) 



It is important to stress that this condition does not 
guaranty the existence of quasi- long-range order at T = 
if gA < igA)crit- To see this one can estimate the 
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number of pairs of sites on which it is energetically fa- 
vorable to create a dipole. This number is approxi- 
mately P(0, +m, —m). In turn, with the choice = o,, 
P(0,xi,X2) can be read off from P(0,xi) provided one 
replaces ln(i/a) in P(0,xi) by 21n(|xi — X2|/a): 



-7rJ(l 



P(0, +m, —to) 
47r 

%tt'^ ,P g Am? 



L 



dy e 



2 1- 



(A17) 



Hence, the number of pairs of sites on which it is en- 
ergetically favorable to create a dipole always diverges 
with system size [like (L/aY if <7a < '^{9A)crit, like 
(L/a)*"^'^"' )/(2sa) otherwise]. Likewise, the number of 
n sites on which it is energetically favorable to create 
a neutral vortex configuration always diverge with sys- 
tem size. Hence, for most realizations of the disorder 
dfj^O, the ground states are not the ferromagnetic state 
= but are neutral and non-trivial [complex) vortex 
configurations. To decide whether such complex ground 
states support quasi-long-range order, one must estab- 
lish absence of screening of the CB potential between to 
external charges. 



2. Screening at zero temperature 

To illustrate the issue of screening, we introduce the 
random energy 

Hgb[Q,Oc] -.^ -T:jJ2(^k - nk)Gkiimi - m), (A18) 

k,l 

where is a neutral configuration of charges to^ — 0, ±1, 
9a is a neutral configuration of charges rij, — 0, ±a, 
< a < 1, that realizes the disorder, and Gki is a 
short hand notation for the logarithmic CB gas poten- 
tial. Finally, we take the probability to realize da to be 
proportional to 



exp 



nkGkini 



(A19) 



The relationship between this model and the random 
phase XY model is that only a very small subset of all 
possible realizations of the disor der in the random phase 
XY model are allowed in Eqs. ( A18 , Al9| ). Thrown out 
are all realizations of the disorder made of vortices of 
unequal vorticity. 

Consider now the vortex configuration defined by 
77ifc = 0, ±1, respectively, whenever = 0, ±a. The en- 
ergy of this configuration, whose unit vortices track the 
fractionally charged quenched vortices, is 



of^mkGkimu 

k,l 



(A20) 



and should be compared to the energy of the ferromag- 
netic state 

i?CB[O,0a] := -^Ja'^TOfcGfeiTO,. (A21) 

k,l 

We see that it is energetically favorable to create a unit 
vortex TO/c at the location of every quenched vortex Uk 
provided 1/2 < a < 1. Otherwise, the ferromagnetic 
state is preferred. For a — 1/2, Sc, and 9a are degener- 
ate ^. 

We are in position to ask the following question. Is 
the CB potential screened or not for the family of vortex 
configurations {Sq} where 1/2 < a < 1 is held fixed? 
Since merely creates vortices wherever quenched vor- 
tices sit, the question can be reduced to: what are the 
screening properties of the CB gas with effective temper- 
ature and charge qa and a, respectively? The answer is 
known p^ ], namely for sufficiently small qa, i.e.. 



9A<a 



(A22) 



the CB gas does not screen since the dipole phase is real- 
ized and there exists quasi-long-range order. Otherwise, 
the CB gas screens since the plasma phase is realized and 
quasi-long-range order is destroyed. 

Once we know the asymptotic form of the CB poten- 
tial between two external charges for the family {Sa}, 
1/2 < a < 1 and gA fixed, we can extract the a-averaged 
CB potential where we restrict 1/2 < a < 1. The a- 
averaged CB potential is controlled in the ther mody - 
namic limit by a = 1/2 since the probability in Eq. ( A19| ) 
scales like exp[— a^(L/a)^/(L/a)], where f{x) is some 
positive function with lima;|oo f{x)/x^ = that does not 
depend on a. The a-averaged CB potential, 1/2 < a < 1, 
thus decays logarithmically for sufficiently small qa and 
decays exponentially if > (3A)crit- Hence, quasi-long- 
range order is present for not too strong disorder strength 
provided it can be shown that Sq, is the ground state for 
every realizations of the disorder when 1/2 < a < 1. 

Although this conjecture might be reasonable for the 
toy model, it is certainly not true for the full problem 
where a given realization of the disorder 8^,9 can create 
vortices of arbitrary fractional charges in contrast to the 
toy model. Hence, although the existence of both 9a and 
Sq, suggests the existence of {gA)cTit in the full model, 
we cannot rule out more complex configurations O with 
energies lower than that of both the ferromagnetic state 
and Sq, and for which the CB potentials are screened at 
long distances for some gA < (gA)crit- 



^We are indebted to E. Fradkin for this observation. 
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APPENDIX B: FUGACITY EXPANSION 

In this section, we we are going to expand all cor- 
relation functions of the SG operator e'^ in powers of 
hi/2t and then perform disorder averaging order by or- 
der in powers of hi/2t. We thus reproduce the fugac- 
ity expansion performed by Korshunov on the CB 
gas with quenched fractionally charged vortices. Scaling 
fields with negative scaling dimensions are then seen to 
control the singular behavior of the expansion in powers 
of hi/2t of the inverse SG partition function order by 
order. However, we begin first by illustrating a possible 
drawback of the fugacity expansion. 



1. Drawback of the fugacity expansion 

The expansion of the inverse SG partition function in 
powers of the fugacity hi/2t should be taken with great 
caution since convergence is not always warranted. In- 
deed, let X be a real valued random variable and let 

Y := 1/(1 -I- X) be another real valued random variable. 
For example, Y could be 1/Zcb whereby X could be 
ZcB — 1- Let Px{x) be the probability that X takes the 
value X. We want to calculate the probability Py{v) that 

Y takes the value y. For definiteness. 



• Case I: 

Px{x) 

• Case II: 

Py{v) = 



e 
1- 



< X < oo 

0<y<l. (Bl) 



We immediately conclude that the expansion 



(B2) 



(B3) 



would predict that all moments of Y diverge in both cases 
I and II. Howeve r, th ese moments can be calculated di- 
rectly from Eqs. (Bl B^ ) and are all finite. 



2. Preview to the fugacity expansion to fourth order 

Let n be a positive integer, choose n points xi, • • •, x„ 
on the Euclidean plane, and define 



F — pieix(xi) . . i£„x(x„) 

^ Xi ,■ ■ - jXti ■ — ^ ^ ; 



jv[x] e'^^°[^^^lFx„ 



(B4) 
(B5) 



Each factor g'^fcxCxt)^ ^ _ . . . ^ can be thought of as 
the insertion of an external vortex of vorticity = ±1 in 
the CB gas. Thermal averaging of Eq. (B4) is obtained 
by dividing the unnormalized average in Eq. (B5) by the 
SG partition function: 



(^x,,.,xj 



(B6) 



Finally, disorder averaging over 9^0 is done with the 
probability distribution of Eq. (3.4). 

We attempt to calculate both (Fxi, ••,x„) and 
(^xi,- -,x„)"""°' through a power expansion in hi/2t and 
then perform term by term disorder averaging over 8^,9. 
Thermal averaging and 9^(0 averaging are seen to "fac- 
torize" to each order in hi/2t for (-Fxi, -.xj""""'- We 
then go on performing d^0 averaging over (Fxi,...,x„)- 

The key identity that is needed is 



(^xi,..,x„)~ = (B7) 

i:;^(|)"'i:C;)/a.-7.v„.x 

m=0 ^ ^ p=0 \ ^ y 

/gix(yi) . . . gix(yp)g-ix(yp+l) . . . g-ix(ym)^_^ ^ \unnor 
\ xi,---,x„/ ^^^^ 

Thermal averaging on the last line must be performed 
with hi — 0, in which case averaging over x is Gaus- 
sian. We can then use the shift of integration variable 
X =: x' + to decouple averaging over x' from aver- 
aging over dfj^O: 



(B8) 



E 



1 



m=0 "" ^ ~ ' p=0 \ ^ / 

g- 2Ve(yi) . . . g2Te(y-)e- 2V ELi ""'^^'"'^ x 

. , , , . , , , . v^ii // ^ \ unnor 

gix (yi) . . . g-ix (y™)gi2^fc^i ^fex i^k) \ 

I /ii=0 



Thus to each order in hi/2t, averaging over d^6 has 
factorized from averaging over x' in the integrand on the 
right hand side of Eq. (B8). In fact since both averages 



are Gaussian, one obtains the CB gas representation 



/ p \ unnor — 

\-'^xi,---,x„/ — 



exp 



p=0 



h 



(B9) 



where 
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Yfe, iffc = l,---,m, 

Xfe, if fc m + 1, • • • , TO + n, 



(BIO) 



and we have introduced the effective CB gas couphng 
constant 



2'kK := 2tt{K - qaK^) = 27r 



47r2 ^ 47r2 I 27r 



(Bll) 



Equations (B9 Bll) tell us that under the assumption 
that a fugacity expansion of the CB gas is vahd for each 
reahzations of the disorder 9^0, then averaging over 5^0 
ah unnormahzed correlation functions consisting in the 
insertion of external charges amounts to the same renor- 
malization of the CB gas effective temperature 1/ttK. 
Moreover, in view of 



E l / IJ^ 
to! V 2t 

m=0 ^ 



, \ m m 

p=0 



(fyi ■ ■ ■ (fy„i X 



e'x(yi) . . . gix(yp)g-ix(yp+i) . . . g-ix(ym) 



(B12) 



this is also true for the case n = which corresponds 
to the average of the partition function expanded in 
powers of the fugacity. We thus expect that the fugac- 
ity expansion for all unnormalized correlation functions 
(^XL - .x,!)"""""^, if well defined, is convergent in the same 
region of the l/K, g^, hi/2t couphng space. However, 
this is not so fo r the no rmalized average (i^xi, --,x„)- 

Indeed Eqs. ( |B9|JbTT1 ) are not sufficient to extract the 
effect of disorder averaging over d^0 on the perturbative 



expansion in powers of hi/2t of (Fx 



In fact, we 



need all moments (not only the first one) of the unnor- 
malized correlation functions (Fxi,....x„)"""" as can be 
seen by expanding 1/2^sg[^] in powers of hi/2t in Eq. 
(pq). But these moments can be estimated from 



:ix(yi) 



=ix(yp)e^'x(yp+i) 



.e-ix(y™)\ — 

^ /hi=0 



exp 



where 



-nK{q) ^ EfcEiln 



k,l=l 



Yk - yi 







Kiq) ■.= Kq-gA{Kqy 



(B13) 



(B14) 



The boundary K{q) = which delimits a positive from 
a negative effective CB gas temperature shrinks with in- 
creasing q. Hence order by order in powers of hi/2t, the 
regime in which the fugacity expansion is well defined is 
controlled b y the largest moment contributing to this or- 
der in Eq. ( B13 ). We recognize in Eq. ( B14 ) the scaling 
dimensions on the right hand side of Eq. (5.7). 



3. CB gas interpretation of the fugacity expansion 

Before going into a more detailed discussion of the fu- 
gacity expansion up to fourth order, we comment on some 
general features of the fugacity expansion and its close 
relationship to the replicated effective theory of section 
0. The impurity average over the two-point correlation 
function (^12) in Eq. ( |4.2| ) can be recast as a summation 
over all possible configurations 



('eix(xi)-ix(x2) X e'X(yi) + ----ix(y2i)\unnor ^ /gix(zi ) + • • •-ix(z23 ) \ unnor ^ 
\ / h-i—0 \ /fii—0 



(B15) 



oi 2i + 2j -\- ■ ■ ■ vortices that screen the two external vortices at Xi and X2. The appropriate combinatorial weight 
that results from expanding the numerator and denominator in powers of the fugacity is not written here but can be 
found in the coming subsection. The first thermal average (• • comes from the expansion of the numerator in 

Eq. (B.15), while the remaining factors in (• • •)JJ""o'^ x • • • come from the expansion of the denominator . Notice that 
thermal averages are disconnected averages as is indicated by the superscript [see Eq. (B5)]. 
Introducing x =■ x' + '^(^ (since x' and decouple), the above equation becomes 



+ ---/gix'(xi)-ix'(x2) X e'X'(yi) + ----ix'(y2.)\unnor x (e^X' + - ■ --^X' {22j))annor ^ 

\ ' lll—U \ ' lll—V 



(B16) 



r 



Each term in the fugacity expansion can then be given 
the following CB gas interpretation. There is one charge 
Q, the disorder charge, associated with the field 9. It 
is also convenient to associate with each disconnected 
thermal average (• • O™"" appearing in ( B16| ) a distinct 
thermal charge £a- More precisely, to order {hi/2t)'^^ 
we introduce n -\- I thermal charges Sa, a = 0, - ■ ■ ,n (we 
take the thermal charges e„i+i = e„ = to al- 

ways vanish if there are only m < n thermal factors, and 



the charges labeled by the subscript always refer to the 
disconnected thermal average involving the two external 
charges at Xi and X2, respe ctively). For example, if none 
of the coordinates in (B16) coincide we assign 



• (eo = +1, £1 = 0, • • • , £„ 0; Q = 

charges at xi and yi, I — 1, - ■ ■ ,i. 

• (eo = -1, £1 = 0, • • • , £„ = 0; Q 

charges at X2 and yi, I = i + 1, ■ 



-1) to the CB 



-1) to the CB 



.,2z. 
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• (eo = 0,ei = +1, •••,£„ = 
charges at zi, I — 1, - ■ ■ ,j. 



0;Q 



-1) to the CB 



• (eo = 0,ei = -1, • • • ,e„ = 0;Q = -1) to the CB 
charges at z;, Z = j + 1, • • • , 2j. 

Since coordinates belonging to distinct disconnected 
thermal averages can coincide (see the next subsection), 
we will also allow complex or fused charges of the form 
(eo,---,e„;(5) where 



^e„ = Q, |ei|<n, |£2| < n - |ei|, • 



(B17) 



a=0 



It is then possible to systematic ally carry through the 
thermal and disorder average in (B16) by assuming that 
the replicated charges (eo, ei, • • • ; Q) and (eq, e'j^, • • • ; Q') 
are associated to operators with correlation functions 
given by 



-2TvgAK'^QQ'+2TvKj2^ ^a^'a 



(B18) 



As we change the microscopic cut-off a and/or in- 
clude randomness in hi, two operators with charges 
(eoi £ij • • ■ ; Q) and (eg, e'j^, • • • ; Q'), respectively, may fuse 
into one with charge {e^ -I- ef,, e^ -I- e'^, • • • ; Q + Q'). Thus, 
operators with arbitrary complex charges may appear in 
the fugacity expansion upon rcnormalization or averag- 
ing over random vortex fugacity. Although bare vortices 
in the XY model have a simple structure (being labeled 
by a single integer), it is striking to see that complex 
vortices must be accounted for in disorder averaged cor- 
relation functions in contrast to the clean limit (in the 
pure system these higher charges are always irrelevant). 
Complex vortices also appear formally in the replica ap- 
proach (see section ^). However, intuition is easily lost 
when taking the replica limit r J, 0. 

Our discussion of the fugacity expansion suggests that 



fugacities Yi 



{eo,ei,---:,Q) 



should be introduced when op- 



erators with charge (eo, ei, • • • ; Q) contribute to ( B16D . 
These fugacities, to a first approximation, are related to 
the density of screening charges of type (eo, ei, • • • ; Q). 
Form invariance of the fugacity expansion under an in- 
finitesimal rescaling of the short distance cutoff a: 



ae 



< Z< 1, 



(B19) 



would then imply that the fugacities obey RG equations. 
Since the scaling dimension of an operator depends on 
its charge, we expect the RG equations for the new fu- 
gacities to be different from each other. In particular, 
to order 2n in the fugacity, the generalized fugacities 



(to order 2n in the fugacity expansion). For example, to 
fourth order, we will denote by y(i,i;2) any of the three 
fugacities y(i,i.o;2)j Y(^i,o,i-2), and y(o.i,i;2) associated to 
the charges (1, 1, 0; 2), (1, 0, 1; 2), and (0, 1, 1; 2), respec- 
tively. Indeed, we will show that the RG equations for 
all three fugacities are the same. 

We illustrate those general considerations with a de- 
tailed calculation of the four th order correction to 

(-Fxi,x2)- We denote with F^2' the fourth order coeffi- 
cient to the fugacity expansion (see the following subsec- 
tion). It is then convenient to distinguish between three 

contributions to , denoted A]^2' ^121 and C12, re- 
spectively. To fourth order in the fugacity expansion, we 
expect that complex vortices with charges (eg, ei, e2; Q) 
emerge. Correspondingly, it should be possible to recast 
the RG equations in terms of fugacities ^(eo.£i.e2;Q)- This 
is indeed so. 

The contribution is nothing but the second order 
contribution squared. Hence, it is due to the screening of 
two external vortices with charges (-f 1; -1-1) and (—1; —1), 
respectively, by four thermal vortices with charges of type 
(±1;±1). It is given by: 



12 I 2t 



X12 



-2wK 



27ra;(2) In 



X12 
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where x^^-* is the second order correction of Eq. (4.5) 
when q — 1: 



L/a 



n(2) 



-4^3 [^2 _ ^4^2 ] y2^ l^yy 



,,3-2ttK 



(B21) 



This is an encouraging result since it justifies a posteriori 
our assumption in Eq. (4.4) that logarithmic corrections 
can be reexponentiated. 

The contribution is due to the screening of our 
two external vortices by either three replicated vortices 
with charges of type (-1-1, -1-1; -1-2) (fused), (—1,0; —1) and 
(0, —1; —1), respectively, or three replicated vortices with 
charges {+1,-1;0) (fused), (-1,0; -1) and (0, -1-1; -1-1), 
respectively. 

The term is induced by the screening of two exter- 
nal vortices at Xi and X2 , respectively, by either two fused 
vortices with charges of type (±1, =f1; 0), or by two fused 
vortices with charges of type (±l,±l;±2). Concretely, 
we find 





X12 




a 



-2-kK 



X 27ra;l'^) In 



X12 



(B22) 



Y, 



(0,1,---,1;Q)' ■ 



4,0;Q) 



are expected to be the most where x^c'' is the fourth order correction due to two-body 



relevant operators within the family y(eo, - >e7i;Q) [com- 
pare with Eq. (5.7)]. 

We will slightly abuse our notation in the following 
by using for the subscript of generalized fugacities the 
maximum number n of non-vanishing thermal charges 



rcnormalization effects: 



L/a 

44):^_47r3if2y(2 _i^g) jdyy 



3~4ttK 



(B23) 
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K{2) 



Y2 



L/a 

1 



The first line on tlie right hand side of Eq. ( B23 ) is 
fully oblivious of the disorder. The second line is dra- 
matically sensitive to it since it yields a new boundary 
for analyticity of the fugacity expansion. The two lines 
are equal in the absence of disorder. 

Fu rther more, the second line on the right hand side of 
Eq. (B23) implies that, in contrast to the second order 



contribution to the fugacity expansion of the mean {F12), 
the scaling exponent 2nK{2) of (-Fi2)^j^o enters in the 
fourth ord er correction to the scaling exponent 2ttK of 
{Pi2)hj^^o- Hence, the second moment of (F12) couples 
to the mean of {F12) in the RG flow to fourth order in 
the fugacity. 

From the scaling dimension 2TrK — Ang^K^ = ttK{2) 
of the operator with charge of type (1,1; 2) , we find 
the RG equation ( [4.16D of Y(i,i;2)- Equation ( p23| ) tells 
us that y(i,i;2) indeed enters the RG equation of K. 
H ence, t he perturbative RG equations breaks down when 
2ttK{2) < 4. 

The pattern for renormalization should become clear 
from our fourth order calculation. To each new order 
2n in hi/2t we need to introduce new fugacities for 
replicated vortices labeled by their thermal and disorder 
charges 



to close the RG equations. The scaling dimension of the 
fugacity 

on the critical plane of Fig. 1 is deduced from that of 



exp 



exp 



a=l 
n 



exp 



(B24) 



For any infinitesimal value of the disorder strength gA , 
the contribution to the coefficient of the fugacity expan- 
sion to or der 2 n that defines the regime of analyticity 
[see Eq. ( 4.14 )] describes the screening of the CB in- 
teraction by two tightly bound fused vortices of charges 
(£,•••,£; ne), e = ±1, respectively. The scaling wit h sys - 
tem size of this coefficient is thus given by Eq. (4.13). 
Subleading contributions to the coefficient of the expan- 
sion originate from screening of the CB interaction by 
three and more replicated vortices. 

Finally, we would like to stress that there exists a one 
to one correspondence between the replicated vortices ap- 
pearing in the fugacity expansion and the replicated pri 



mary fields constructed in section ^ see Eqs. (5.5,^^)] 



as was first suggested by Korshunov [ 2ll 



(eo,- 



D. Two-point function (Fxi,x2) up to order (hi/2t) 



It is very instructive to carry the fugacity expansion of the two-point function 

P — pix(xi)-ix(x2) — p,„ 



(B25) 



This calculation is the crucial ingredient in the RG analysis of the KT transition. 

For the pure case quasi-long-range order holds if F12 decays algebraically with separation with a scaling exponent 
K = x{l) [see Eq. ( |4.4| )] which is an analytic function of fugacity Yi in the vicinity of Yi — 0. The transition 
temperature to a disordered phase is deduced from the boundary along the line 1/K > for which the scaling 
exponent k becomes non-analytic. 

In the presence of disorder, Rubinstein et al. performed the same analysis after averaging over disorder the fugacity 
expansion of F12 term by term up to second order in Yi. They inferred the parabolic boundary of Fig. 1 from 
the onset of a non-analytic dependency of k on fugacity Yi. We repeat their argument and show, in the spirit of 
Korshunov's analysis [|21[| , how higher moments of correlation functions invalidate their conclusion beyond second 
order in the fugacity expansion. 

The power expansion in Yi cx hi/2t is given by 



12} 



n=0 



12 



(B26) 



where i^^f = fg^+'^ 



= = and 



F, 



(0) 



12 



p(2) 
^12 



F, 



(4) 



12 



_ AO) 

— J12 ' 

_ A2) _ AO) y{2 

— J12 Jl2 ^ 

_ f(4) 
~ il2 



(B27) 
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up to fourth order in hi/2t. The coefficients in the power expansions in hi/2t of the numerator and denominator are 



„(2n) 
J 12 



1 



(n!)2 
1 



d^Yi ■ ■■d^y2n (e 



,ix(yi)H hix(yTi)-ix(yn + i) ix(y2,0+ix(xi)-ix(x2) 



(fyi ■ ■ ■ (fy2n ( e 



,ix(yi)H — i-ix(y.i)-ix(yn+i) ix(y2n) 



(B28) 
(B29) 



respectively. Since we are assuming the existence of a dipole phase, thermal vortices are taken with a hardcore as 
is indicated by the constraint that the coordinates of the vortices cannot coincide. Hence, we are implicitly using a 
short distance cutoff a for the thermal vortices. As a matter of principle, this cutoff need not be the same as that used 
at short distances for the quenched vortices. Nevertheless, for notational simphcity, we will assume that quenched 
vortices share the same hardcore radius. 
To lowest order in Yi c>c hi/2t: 



5^(0) 
12 





-2tvK 




Xi - X2 




Xl2 


a 




a 



'2ttK 
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To second order in Yi cx hi/2t: 



(2)_ 



a 


2 


2itK 


xi - 


- X2 






2-itK 




Xl2 




J^12^ 









d yid y2 



2 ^xix2(yi,y2;27rX) -/Cxix2(yi,y2;27ri4'2^A) 



|yi - y2 



\2ttK 



^xiX2 (yi , y2 ; 2t:K) ~ /Cxix^ (yi , y2 ; 2'KK'^gA) 



|xi2/a| 



|yi2 



\2ivK 



where 



^xiX2 (yi,y2;a;) 



d|yi2| |yi2p|yi2/c 



|yi -xi||y2 -X2I 



X In |xi2/a| 



|yi -X2||y2 -xi| 



(B31) 



(B32) 



If we interpret the presence of a logarithmic correction on the right hand side of Eq. (B31) as the first term in the 
expansion of 



F12 



X12 



X12 



2TrK+SK. 



in powers of Sk, then the scaling exponent k governing the algebraic decay of F12 is given by 



dx X' 



,3-27rA' 



(a4ri2)+0(a8lf). 



(B33) 



(B34) 



We see that the condition that the scaling exponent for the algebraic decay of the two-point function be analytic in 
the fugacity Y\ in the vicinity of vanishing fugacity yields the parabolic boundary in Fig. 1. 

To fourth order in Y\ cx hi/2t we must distinguish between three different effects present both with or without 
disorder. Indeed, we can write 



f[$ = A12 + B12 + C12 



(B35) 



Renormalizations of the interaction between two external vortices due to: i) four-body effects, ii) three-body effects, iii) 
two-body effects, are denoted by A12, B12, and C12, respectively. Three-body and two-body effects must be accounted 
for when the separation |xi2| between the two external charges is much larger than the hardcore radius a, and when 
the separation between vortices is within the hardcore radius. Formally, we implement these renormalization effects 
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whenever coordinates of vortices coincide in the integrands. This is never allowed due to the hardcore constraint 

f or b oth /^Ij"* ^^'^ friZ^^h However, coordinates can coincide when expanding the inverse partition function in Eq. 
( B26D . For example, 



Al^Zi^)= J £yi£y3 J £j^£j3(eUx(y.)-x(y3)+x(x.)-x(x.)])— x (e'[x(y2)-x(y4)]) — 

5^ ^ 



1234 



/ei[x(yi)-x(y3)+x(xi)-x(x2)]\"""°'^ x /ei[x(y2)-x(y4)]\™" 

\ / hi—0 \ I /ii— 



+ / (ei[x(yi)-x(y3)+x(xi)-x(x2)])— x (eiMyi)-x(y4)]) — 
+ / /ei[x(yi)-x(y3)+x(xi)-x(x2)]\"""°'' x /e'[x(y2)-x(yi)]\"""°' 

I \ I hi—O \ I hi—0 

Jl32 

+ / /e'[x(yi)-x(y3)+x(xi)-x(x2)]\"""°'' x /e'[x(y3)-x(y4)]\"""°' 

I \ I /ii— \ / /ii— 

J 134 

+ I /ei[x(yi)-x(y3)+x(xi)-x(x2)]\™"" x /ei[x(y2)-x(y3)]\"""°'' 

I \ /hi —0 \ /hi —0 

Jl32 

+ f /ei[x(yi)-x(y3)+x(xi)-x(x2)]\™" x /e'[x(yi)-x(y3)]\"""°'' 

I \ I /ii— \ / /ii— 

13 



+ / /ei[x(yi)-x(y3)+x(xi)-x(x2)]\™" x /e'[x(y3)-x(yi)]\™°'' 

f \ I hi—0 \ I /ii— 



(B36) 



Four-body renormalization effects are 
1 



a ^ a |yi2||y34| ^ |yi - xi I \y2 - xi I |y3 - X2 1 |y4 - X2 1 
(2!)^ A234 L 1^12 1 Iyi3||yi4||y23||y24| |yi -X2||y2 -xsHys -xi||y4 -xil 



2ttK 



1 



(2!)2 



3412 



a a"|yi2||y34| 



|xi2| |yi3||yi4||y23||y24| 



27rA' r 



|yi - xi I |y2 - xi I |y3 - X2 1 |y4 - X2 1 
|yi - X2 1 |y2 - X2 1 |y3 - xi I |y4 - xi I 



3412 
1432 



a a 



|xi2| |yi3||y24| 



2ttK 



|yi2||y34| |yi -XiHya -XiHyg -X2||y4 -X2I 



X 



|yi4||y32| |yi - X2||y2 - X2||y3 - Xi||y4 - Xil 



a a |yi - xi||y3 - X2I 



-,2ttK r 



|Xl2| |yi3||y24| |yi - X2||y3 - Xi| 



|yi2||y34| |y2 - xi||y4 - X2I 
|yi4||y32| |y2 - X2||y4 - xi| 



h27r/4'^ 



9a 
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Assuming screening by a dilute gas of dipoles reduces the integration over the coordinates yi, • • •, y4 of the four 
thermal vortices to 
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23 



^(13) 
^(21) 
^(23) 



^3 



yi2|>a |yi3|>a 



yi2 



27riC(2)/2 



yi3 - yi2 



a 



yi3 



2irK(2)/2 



yi2 • yi3 



(Pyi3 



J a" 

|yi2l>a |yi3|>a 



(Pyi2 f (Pyi3 



yi2 



2irK 



yi3 - yi2 



a 



+27rK^g^ 



yi3 



2irK 



yi2 • yi2 



a" J 

|yi2|>a |yi3|>a 



yi2 



2irK 



yi3 - yi2 



a 



+2TrK^g^ 



yi3 



27tK 



yi2 • yi3 
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Here, we have introduced the dimensionless fugacities y(ei,e2;£i+£2)' ^(o,e;£)j ^^nd Y(^^q.^-j, where £i = £2 = il and 

/ 2, \4/3 

e = ±1, respectively. Their bare values are given by ( ^2t ) • Whereas under a RG rescaling a = e'a' with 
< Z <C 1, the y's and product thereof are renormalized niultiplicatively, an open problem is to find numerical factors 
entering additive renormalization effect due to three-body effects. This is difficult in this real space RG approach 
since we must perform four-dimensional integrals with complicated integrands and intricate boundaries. 

Form invariance under an infinitesimal rescaling of the short distance cutoff a: a' := ae', < ^ 1, of B12 implies 
that the three fugacities {ei,e2 — ±1) 



(B48) 



obey: 



1. They always appear in the combinations 

^(£i,£2;£i+e2) ^ ^(-£i,0;-£i) X "^(0,-£2;-e2)- 



(B49) 



2. Those combinations rcnormalizc multiplicatively according to two rules (e ~ ±f ) 

[5^(s.e;2.)y(-.A-e)^(0,-e;-e)]' = l^(e,.;2e)^(-e,0;-s)>^(0,-s;-e) e6-2-^-2.K's 



Additionally, if we require the individual multiplicative renormalization rules 



Y' 



(e,£;2e) 



Y, 



J[2-kK(2)]^ 
J[2-2ttK\ 



Y! 



(£,0;e) 



^Y, 



Y' 



(0,£;£) 



= Y, 



Jl2-7rK)]^ 
J12-7TK] 



we can summarize Eqs. ( B50 , B50|) by 



^(+l, + l:2) X ^(-1,0;-1) 



X Y, 



(0,-l;-l) 



Y, 



(+1 -1;0) 



X >?-1.0:-l) X Y 



(0, + l; + l) 



gi(xi+X2)(yi) y, g-ixi(y2) y, g-ix2(y3) 
gi(xi-X2)(yi) y, g-ixi(y2) y, g+ix2(y3) 



(B50) 
(B51) 
(B52) 
(B53) 

(B54) 
(B55) 



Equations (B54,B55) tell us that any of the three thermal vortices renormalizing in B12 the CB interaction between 
two external vortices can be thought of as some local linear combination of the two replicated real scalar fields Xi,2- 
We remember that xi.2 couple to the disorder in such a way that the fields x'l 2 defined by Xi 2 — xi 2 + i^ttKO are 
free scalar fields independent of the disorder 9. The fugacity Yf^^^^^^.^^+^^-j is thus labeled by three charges (measured 
in the appropriate units): 

• The thermal charge ei of Xi- 

• The thermal charge £2 of X2- 

• The disorder charge £1 + £2 of 9. 
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Finally, two-body renormalization effects are 
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a 




i|xi2U 



2TrK 



(B56) 
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2x2tiK 
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|yi -xi||y3 -X2I 
|yi -X2||y3 -xil 



2TrK 



The major difference between C12 and Eq. (B31) is the exponent of the coordinate yi3 given by K{2) ~ 2K — AK^ g a. 
It is possible to show that 



^12 







2ttK 
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In 


X12 
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a 
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Here, we have introduced the dimensionless fugacities ^(i.i;2) and y(i,-i;o) whose bare values are equal to (^^2!^^ ■ 

In the absence of disorder, renormalization effects due to higher charge vortices do not modify the boundary 
extracted from Eq. (B31). This is nothing but the statement that cos(x) is the most relevant operator of the family 
cos(xi) X • • • X cos(xg), • • • ,cos(gx), q G N, along the Gaussian fixed hue 1/K > 0, gA — 0, Yi — 0. However, in the 
presence of the random vector potential d^9 the relevance of cos(xi) x • • • x cos(xg) increases with q e N. 
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